Hydrodynamics of Bose and Fermi superfluids at zero temperature: 
the superfluid nonlinear Schrodinger equation 
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OO ' We discuss the zero-temperature hydrodynamics equations of bosonic and fermionic superfluids 

^^ \ and their connection with generahzed Gross-Pitaevskii and Ginzburg-Landau equations through a 

>— ' . single superfluid nonhnear Schrodinger equation. 

CN . 

^ ' PACS numbers: 03.75.Lm, 03.75.Ss, 05.30.Jp, 74.50.-|-r 

o' 

^H ' Recent and less-recent experiments with ultracold and dilute gases made of alkali-metal atoms have clearly shown 

the existence of superfluid properties in these systems 1, 2]. Both bosonic and fermionic superfluids can be accurately 
described by the hydrodynamics equations of superfluids [ll, 0, Q- I^^ this paper we analyze the hydrodynamics 
equations of superfluids and show how to construct a reliable nonlinear Schrodinger equation from these hydrodynamics 
' (-H ', equations. For bosons this equation gives the generalized Gross-Pitaevskii equation recently discussed by Volovik [4|, 
^ ' while for fermions one gets a zero-temperature Ginzburg-Landau equation [5] . The limits of validity of these mean-fleld 
O , equations are discussed. 

At zero temperature the hydrodynamics equations of superfluids made of atoms of mass m are given by 
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where n(r, t) is the local density and v(r, t) is the local superfluid velocity [ll, l2, Q- Here U{y) is the external potential 
■ O ■ and /i(n) is the bulk chemical potential of the system. The bulk chemical potential /i(?T.) is the zero-temperature 
C , equation of state of the uniform system. The density n(r, t) is such that 



N = f n{r,t) d^r (3) 



^ ' is the total number of atoms in the fluid. Eq. ([IJ and ([2]) are nothing else than the Euler equations of an inviscid 

t — I (i.e. not-viscous) and irrotational fluid. In fact, at zero temperature, due to the absence of the normal component, 

C ■ 1 the superfluid density coincides with the total density and the superfluid current with the total current [3|. 

^^ ' The condition of irrotationality 

t:^' i V A V ^ (4) 

o ; 

QQ I means that the velocity v can be written as the gradient of a scalar field. Eqs. ([T|) and (O differ from the corresponding 

^^ ■ equations holding in the coUisional regime of a non superfluid system because of the irrotationality constraint ((3]). In 

*~J ' addition, experiments with both bosonic and fermionic superfluids show the existence of quantized vortices, such that 

• i-H the circulation C of the superfluid velocity 

^ . C=fx-dv, (5) 



is quantized, i.e. 



C = —: , (6) 



where k is an integer quantum number and the statistical coefficient C is 1 for superfluid bosons and 2 for superfluid 
fermions [H, 0, i, [ll . 

Eq. ([6]) does not have a classical analog, and this fact suggests that the superfluid velocity is the gradient of the 
phase 0(r, t) of a single- valued quantum-mechanical wave function S(r, t). This function 

S(r,t) = |S(r,t)|e'^(-->*) (7) 



is the so-called macroscopic wave function of the Bose-Einstein condensate of the superfluid [3, [I4I . The connection 
between superfluid hydrodynamics and quantum field theory is made by the formula 

V = Ave , (8) 

Qm 

where again C = 1 foi' bosons and C, — 2 for fermions. For bosons the condensate wave function is given by 

S(r,<) = (^(r,t)), (9) 

that is the thermal (or ground-state) average of the bosonic field operator ipir^t) [9|]. For fermions the condensate 
wave function is instead 

E{v,t) = {■i^{v,t)i,y{v,t)) , (10) 

that is the average of pair operators, with '(/'^(r, t) the fermionic field operator with spin component a =t , | [9|. Notice 
that the condensate wave function S(r, i) is normalized to the total number of condensed bosons [9| or of condensed 
Cooper pairs [l2, [l3| -, and in general the condensed fraction can be much smaller than one [^, ll^ • 

The effect of statistics enters in the explicit formula [ll, SSl of the bulk chemical potential /z(n) which appears in 
Eqs. ([T]) and ([2]). For instance, in the case of a dilute and ultracold wekly-interacting Bose gas one has 

, , 47rft2 / 32 ^ , \ 

^i{n) = n 1 + —^aBU^'^ , 11 

where as is the inter-atomic Bose-Bose s-wave scattering length and asn^'^ is the gas parameter of the bosonic Lee- 
Yang expansion [7| . Very recently we have proposed a Pade approximant to describe iJ,{n) from the the weak-coupling 
regime, where asn^^^ <C 1 and Eq. (fTTj) holds, to the unitarity limit, where a -^ +00 [10]. In the case of a dilute and 
ultracold two-component weakly-interacting Fermi gas one can take instead 

^("^ ^ £ (^"'"^'^' (^ + ^ ^^^'^"' "^"'^') ' (^^^ 

where ap is the inter-atomic scattering length of atoms with different spin and apv}'^ is the gas parameter of the 
fermionic Huang- Yang expansion |8|. Few years ago we proposed a reliable fitting formula of /i(n) based on Monte 
Carlo data in the full BCS-BEC crossover from weakly-interacting Cooper pairs, where Eq. P^ holds, to the Bose- 
Einstein condensate of molecules [ll| . 

The hydrodynamics equations ^ and ([2]) are valid to describe the long-wavelength and low-energy macroscopic 
properties of both bosons and fermions. In particular, one can introduce a healing (or coherence) length S, such that 
the phenomena under investigation must be characterized by a wave length A much larger than the healing length, 
i.e. 

A>e. (13) 

As suggested by Combescot, Kagan and Stringari [15], the healing length can be defined as 

^=^, (14) 

mVcr 

where Vcr is the Landau critical velocity above which the system gives rise to energy dissipation. For bosons the 
critical velocity coincides with the first sound velocity, i.e. 



-- = \1-'^- (15) 

V m on 

For fermions the critical velocity is instead related to the breaking of Cooper pairs through the formula 



"=^ = V m ' ('^^ 

where JA] is the energy gap of Cooper pairs [3, Il3- We notice that in the deep BCS regime of weakly interacting at- 
tractive Fermi atoms (corresponding to ] A] ^ /i) Eq. p6p approaches the exponentially small value Vcr = ] Aj/^/SmTt. 



In addition, we remind that some years ago Kemoklidze and Pitaevskii derived the zero-temperature Eqs. ([T]) and ([2]) 
for the BCS Fermi gas starting from the Gorkov equations of quantum-field theory under the assumption of neglect- 
ing the spatial derivatives of the energy gap |A| [1^. In the case of a Fermi gas that performs the whole BCS-BEC 
crossover the critical velocity Vcr can be estimated as the minimum value between the two values given by Eq. (jlSp 
andEq. ^ [15,]. 

Inspired by the Ginzburg-Landau theory of superconductors [l7| and by the density functional approach to the 
superfluid ■*He [l3l, we now try to express Eqs. H]) and ^ in terms of a Schrodinger equation such that its wave 
function '^{r,t) is given by 



,I,(r,i) = ,/!^Me^«(M) 



(17) 



where the coefhcient C is 1 for bosons and 2 for fermions. We call C the statistical coefficient. In this way the function 
^(r,i) describes superfluid bosons (C = 1) or boson-like Cooper pairs {( = 2) with the normalization 



J\'i{r,t)\^dh^j, 



(18) 



that is quite different from the normalization of the condensate wave function S(r, t) [3, [T^l- Nevertheless, the phase 
9{r, t) of the the complex field ^(r, t) is the same of S(r, t) and also of the gap function A(r, t). Obviously this phase 
must satisfy Eq. ([8]). 

We consider the following nonlinear Schrodinger equation 
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where a, [3 and 7 are parameters which must be determined. By inserting Eq. (|17p into Eq. (|19p . after some 
calculations and taking into account Eq. ([H]), we find two hydrodynamics equations 
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which include the quantum pressure term 
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which depends explicitly on the reduced Planck constant h. This term is necessary in a realistic superfiuid model 
to avoid unphysical phenomena like the formation of wave front singularities in the dynamics of shock waves [l9|. 
Neglecting the quantum pressure term (which is small for a large number of particles apart very close the surface) , 
one gets the classical hydrodynamics equations, i.e. Eqs. ([T]) and ([2]), only by setting 



1 
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/3 = 7 = C 



(23) 



In conclusion Eq. p^ becomes 
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that we call superfluid nonlinear Schrodinger (SNLS) equation, and the quantum pressure term reads 
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For superfluid bosons C — 1 ^^^d from Eq. (f24|) we get a generalized Gross-Pitaevskii equation [J, [lO, l20| which 
becomes the familiar Gross-Pitaevskii equation if /z(n) = (47r?i^as/»Ti)n [Ij|. For superfluid fermions one has instead 
^ = 2 and from Eq. ([M|) we get a zero-temperature generalized Ginzburg-Landau equation [^, [2l| . We stress that 



the appearance of the statistical coefficient ( in Eq. ([24]) is a direct consequence of the relationship ([8]) between the 
phase and the superfluid velocity, while the normalization ([T51) does not affect Eq. ([M)) . 

It is important to stress that Eq. (fT9|) is not the more general Galilei-invariant Schorodinger equation. In fact, we 
have recently shown [23| that a nonlinear Schrodinger equation of the Guerra-Pusterla type |24| , 
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is needed to accurately describe the surface effects of a ultracold superfluid Fermi gas with infinite scattering length 
in a harmonic trap [23j . 

Eq. (|24p is the Euler-Lagrange equation of the following Lagrangian density 
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is the bulk energy per particle of the bosonic superfluid, namely 
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Our SNLS equation (|24p can be used to study stationary configurations and elementary excitations of the superfluids. 
To obtain the stationary equation we set 



*(r,t) = *eg(r) e 



iC,^t/h 



where \I'eg(r) is the equilibrium wave function which satisfies the stationary equation 
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with neg(r) = Cl^eg(r)P the stationary density profile of superfluid bosons {(, — 1) or fermions {Q = 2), and p. the 
chemical potential of the inhomogeneous system. Neglecting the gradient term (Thomas-Fermi approximation) from 
the stationary SNLS equation we obtain the same algebric equation for both bosons and fermions, namely 



U{T)+^lineq{r)) =/i, 

from which we get the stationary density profile 

neq{r) ^ fi^^{fl-U{r)) , 



(32) 
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where fi~^{y) is the inverse function of /i(n). This is exactly the equation one finds in the stationary case fronr Eq. 
(PI). Obviously this stationary density profile strongly depends on the shape of /i(n), i.e. it depends on statistics and 
interaction strength. 

We can also study small deviations from the equilibrum configuration *ei3(r) by setting 



*(r,t) = (*,,(r)-K(/.(r,t)) e-<^*/^ 
It is strightforward to show that, after introducing the sound velocity of the bulk system 
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from the linearization of Eq. (|24p the perturbation 0(r, t) satisfies the equation 
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In the uniform case, where U{r) — and so fi — fj.{neq), we set 

0(r, t) = ^e»(k-r-c:^t) ^_ 5g-^(k.r-t^t) ^ (37^ 

and from Eq. ([36]) we find the Bogoliubov dispersion relation 



^ = J ^5;^ U;^ + 2Cmc ne, 2 . 38 

In the long wavelength limit we find the phonon-like spectrum 

to — c{n^q) k (39) 

of sound waves, while in the short wavelength limit we obtain the particle-like spectrum 

2(m 

that is the kinetic energy of a boson if C = 1 or that of a Cooper-pair if ^ = 2. Clearly Eq. ([3D]) is not reliable because 
the superfluid equations are valid only in the low-energy and long wavelength regime. Notice that by considering the 
next-to-leading term in the small-momentum expansion of Eq. ([38]) one gets 

" = ^("-^)^+ 8Cw1k,) ^'' (4^) 

and this phonon dispersion relation distinguishes between superfluid bosons (C = 1) and superfluid fermions (C — 2) 
also if the sound velocity c{neq) would be the same. 

In conclusion we stress that the superfluid nonlinear Schrodinger (SNLS) equation we have introduced can be used 
for both superfluid bosons and fermions in all the situations where the characteristic wavelengths of the phenomenon 
under investigation is larger than the he aling length given by Eq. (|14p. The SNLS equation can be applied to 
study not only collective modes [ij, [2j llll 123 and free expansion [ll, 12j ^M ' but also quantized vortices where the 
characteristic length is the vortex-core size [l|, [lOJ , and tunneling phenomena where the characteristic length is the 
tunneling penetration depth [5|, [27[ . In addition, with the SNLS eq uation one can investigate interesting nonlinear 
effects, like solitons i25i, [2^, shock waves [l^, and also chaos [20|. Finally, we observe that the SNLS equation 
satisfies the requirements suggested by Greiter, Wilczek and Witten [30] to have a well-founded theory of neutral 
superconductors: it is Galilei invariant, it mantains the current-momentum algebric identity, and at low-energy it 
shows a Nanibu-Goldstone boson field with linear dispersion relation, i.e. the phonon spectrum of Eq. p9p . 
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